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The resistance R vs perpendicular external magnetic field H was measured for superconducting 
Nb thin-film microbridges with and without microholes [antidots (ADs)]. Well below the transition 
temperature, integral R{II) measurements of the resistive transition to the normal state on the 
plain bridge show two distinct regions, which can be identified as bulk and edge superconductivity, 
respectively. The latter case appears when bulk superconductivity becomes suppressed at the upper 
critical field Hc2 and below the critical field of edge superconductivity Hc2, ~ 1.7 Hc2- The presence 
of additional edges in the AD bridge leads to a different shape of the R{II) curves. We used 
low-temperature scanning laser microscopy (LTSLM) to visualize the current distribution in the 
plain and AD bridge upon sweeping H. While the plain bridge shows a dominant LTSLM signal 
at its edges for H > ifc2 the AD bridge also gives a signal from the inner parts of the bridge 
due to the additional edge states around the ADs. LTSLM reveals an asymmetry in the current 
distribution between left and right edges, which confirms theoretical predictions. Furthermore, the 
experimental results are in good agreement with our numerical simulations (based on the time- 
dependent Ginzburg-Landau model) yielding the spatial distribution of the order parameter and 
current density for different bias currents and H values. 

PACS numbers: 74.25.F-, 74.25.0p, 74.25.Dw 



I. INTRODUCTION 

The concept of localized superconductivity in bulk 
superconductors was introduced in 1963 by Saint- 
James and de Gennes^. They demonstrated that 
superconductivity in a semi-infinite sample with an 
ideal flat surface in the presence of an external mag- 
netic field H (with amplitude H) parallel to its sur- 
face can survive in a thin surface layer, even above 
the upper critical field i^c2, when bulk superconduc- 
tivity is completely suppressed. Based on the phe- 
nomenological Ginzburg-Landau theory, the critical 
field Hc3 for surface superconductivity, localized near 
superconductor /vacuum or superconductor /insulator 
interfaces, can be calculated as^i^ 



1.695 H,2 ^ H,3 = Hi°^ (1 - T/Teo) 



(1) 



where H^^^ is the upper critical field for surface su- 
perconductivity at temperature T = 0, and Tco is the 
superconducting critical temperature for = 0. This 
theory predicts that in the regime of the surface su- 
perconductivity the order parameter wave function ^ 
decays exponentially with increasing distance from the 
surface on the length scale of the coherence length ^. 

Experimental evidence for surface superconduc- 
tivity has been found by dc transport^— or in- 
ductive measurement^ shortly after the theoreti- 
cal prediction^ Later on, other methods such as 



ac-susceptibility and permeability measurements^"— 
magnetization measurements^*^ surface impedance 
measurement^ and tunneling spectroscopy^^ con- 
firmed the existence of surface superconductivity 
when H was applied parallel to the surface. The 
evolution of the resistance R vs depending on 
the orientation of H relative to the surface was also 
investigated.^ While two different regions for bulk and 
surface superconductivity were clearly observed for 
fields parallel to the surface, no signature for sur- 
face superconductivity was observed when H was 
applied perpendicularly. The in-plane-field depen- 
dence of the critical current Ic{H) in the regime of 
surface superconductivity for H parallel to the bias 
current was described by Abrikosov^^ and studied 
experimentally J^— Park described theoretically the 
evolution of Ic{H) in the state of surface supercon- 
ductivity when the in-plane field H is applied per- 
pendicular to the bias current flow^^. He predicted 
an asymmetry in the critical surface current, result- 
ing from the superposition of surface screening cur- 
rents and external currents. Such an asymmetry has 
not been observed experimentally yet. 

Similar to surface superconductivity, localized su- 
perconductivity can also nucleate near the sample 
edge in a thin semi-infinite superconducting film, 
in a thin superconducting disk of very large diam- 
eter or around holes in a perpendicular magnetic 
fieldj^ — It should be mentioned, that surface super- 
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conductivity and localized states at the sample edges 
in perpendicular field [called edge superconductivity 
(ES)] are qualitatively and quantitatively the same. 
While surface superconductivity has been investigated 
in several compounds like Pb-based alloys^A^i^, Nb 
and Nb-based alloys^^^ polycrystalline MgB2^, 
Pb^ii^i^, UPts whiskers^, NbSe2^\ experimental 
studies on ES in thin film structures are rare^i^. Re- 
cently, the first real space observation of ES was ob- 
tained by scanning tunneling microscopy on Pb thin 
film islands^. 

Localized states do not only occur at sample bound- 
aries but can also be induced by an inhomogeneous 
magnetic field as it appears e.g. above domain walls in 
superconductor /ferromagnet hybrids. This localized 
state is therefore called domain wall superconductivity 
(DWS)^ — . Recently, a Pb/BaFei20i9 superconduc- 
tor / ferromagnet hybrid has been investigated by low- 
temperature scanning laser microscopy (LTSLM) and 
the inhomogeneous current distribution of the sample 
in the DWS state has been visualized^^ LTSLM is 
therefore a valuable tool to visualize the redistribu- 
tion of the current in the crossover from bulk to edge 
superconductivity. 

In this paper we present our investigations on the 
evolution of edge superconductivity in plain and anti- 
dot Nb microbridges in perpendicular magnetic field. 
Measurements of R{H) were performed to compose 
an experimental phase diagram and to identify the 
regions of bulk and edge superconductivity. Then we 
use LTSLM to visualize the current distribution at 
the transition from the superconducting to the nor- 
mal state in both bridges. In addition, we used a 
time-dependent Ginzburg-Landau model to compare 
our experimental findings with theoretical predictions. 



II. SAMPLE FABRICATION AND 
EXPERIMENTAL DETAILS 

A Nb thin film with thickness d = 60 nm was de- 
posited on a single crystal AI2O3 substrate (r-cut 
sapphire) at T ^ 800° C using magnetron sputter- 
ing. Two Nb bridges with width W = 40 /im and 
length L = 660 /im were patterned by e-beam lithog- 
raphy and reactive ion etching into a bridge geometry 
as shown in Fig. [TJ One of these bridges was pat- 
terned with circular microholes [antidots (ADs) with 
580 nm diameter] in a triangular lattice with a period 
of 1.5 /im. 

The samples were electrically characterized in a He- 
lium cryostat at 4.2 K< T <10K and \H\ <20kOe us- 
ing a conventional four-terminal scheme (cf. Fig. [1]). 
For both investigated Nb microbridges we found 
Tco =8.5 K. H was always applied along the z- 
direction, i.e. H = Hcz was perpendicular to the 
thin film surface and the applied bias current /. We 
performed isothermal measurements of voltage V{I) 



characteristics for different T and H values out of 
which we determined the dependence of the dc re- 
sistance R = V/I on H. The data presented be- 
low were obtained with /=lmA, unless stated oth- 
erwise. This corresponds to a bias current density 
J = I/dW ^40kA/cm^. 




FIG. 1: Optical image of the plain 40 /im wide Nb bridge. 
The contact pads used for / and V are indicated. 

To visualize the current distribution for different 
bias points in the T — H phase diagram, we used 
LTSLM.^^^ ii For imaging by LTSLM, the sample was 
mounted on a cold finger of a Helium fiow cryostat, 
which is equipped with an optical window to enable ir- 
radiation of the sample in the (x, y) plane by a focused 
laser beam with beam spot diameter ^ 1.5 — 2 /um^^^ 
The amplitude modulated laser beam (at frequency 
/ ^ 10 kHz) induces a local increase of temperature 
centered at the beam spot position (xq, ^0) in the sam- 
ple. During imaging, the Nb bridge is biased at a 
constant /, and the beam-induced change of voltage 
Ay(xo, ^0) is recorded by lock-in technique as a func- 
tion of the beam coordinates (xo,yo)- The LTSLM 
voltage signal can be interpreted as follows: If the ir- 
radiated part of the sample was in the normal state 
with resistivity pn^ the laser beam induces a very small 
voltage signal AV oc dpn/dT. However, if the irradi- 
ated part of the bridge took part in the transfer of a 
substantial part of the superconducting currents, the 
beam-induced suppression of superconductivity might 
switch the whole sample from a low-resistive state to 
a high-resistive state. Details of the LTSLM signal 
interpretation can be found in Refs. (sslfilj. 

III. RESULTS AND DISCUSSION 

A. Magnet oresistance data and 
Ginzburg-Landau simulations 

Figure [2] shows R{H) measurements of the resistive 
transition at T=4.2K for different values of / for the 
plain [Fig. EJa)] and the AD bridge [Fig. Efb)]. Ah 
curves are normalized to the normal state resistance 
Rn at i^=9.0kOe. Except for the AD bridge at the 
highest current value of 10 mA, all R{H) curves reach 
Rn at the same field value |i^| ^ 8kOe. However, we 
observe a pronounced dependence of the shape of the 
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R{H) curves on /, which we describe and discuss in 
the following. 

For the plain bridge [cf. Fig. [2fa)], at the highest 
current value of 10 mA, we observe with increasing 
|i^| an onset of dissipation (appearance of a finite R) 
at ^3kOe. Upon further increasing |i^|, the slope 
dR/d\H\ steadily increases, yielding a rather steep 
R{\H\) transition curve up to ^O.dRn- At ^0.9Rn 
{\H\ ~4.4kOe) a kink in R{\H\) appears, i.e. with fur- 
ther increasing the slope dR/d\H\ is significantly 
reduced. Upon reducing /, the field value where the 
kink appears stays almost constant; however, the re- 
sistance at the kink steadily decreases, and becomes 
zero for / <0.5mA, i.e. the kink disappears. Similar 
shapes of the R{H) curves (including the above de- 
scribed kink) and their current dependence as shown 
in Fig. [2fa) for the plain bridge have been found in 
when H was applied parallel to the sample sur- 
face. 

The R{H) measurements of the AD bridge 
[cf. Fig. [2jb)] show similar behavior upon variation 
of / as compared to the plain bridge in the following 
sense: Within the same range of (high) bias currents, 
the onset of dissipation appears almost at the same H 
value (for the same value of /) as for the plain bridge. 
Upon further increasing a similar steep transition 
(with slightly smaller slope as for the plain bridge) ap- 
pears, up to the kink in R{\H\)^ which is also present 
for the AD bridge within the same range of (high) bias 
currents. 

However, we also observe distinct differences by 
comparing the AD and plain bridge: The resistance at 
the kink is lower for the AD bridge. This deviation in- 
creases with increasing bias current. Furthermore, for 
the two highest / values the AD bridge shows a sec- 
ond kink in the R{H) curve where the slope dR/d\H\ 
suddenly increases again with increasing H; this fea- 
ture is absent at smaller / and is not seen for the plain 
bridge for all values of /. Finally, for the two lowest / 
values, the onset of dissipation (upon increasing \H\) 
is shifted to larger \H\ values for the AD bridge, as 
compared to the plain bridge. 

In the following, we present an interpretation of the 
R{H) curves described above, starting with the dis- 
cussion of the results obtained for the plain bridge. 
At the highest 7=10 mA, upon increasing the external 
magnetic field from i7 = 0, vortices will enter the sam- 
ple when H is larger than the field of first vortex entry, 
which is rather small for thin-film structures in per- 
pendicular magnetic field. The onset of energy dissi- 
pation can then be attributed to the onset of motion of 
vortices, when the bias current density J exceeds the 
depinning current density Jdpin at a given T and H. 
In this case, upon further increasing i7, the fiux fiow 
resistance will strongly increase, i.e. , the rather large 
slope dR/dH should correspond to the bias-current- 
stimulated motion of the vortex lattice in the presence 
of a strong pinning potential. The kink in the R{H) 




H (kOe) 

FIG. 2: (Color online) R(H) curves (normalized to normal 
state resistance Rn) measured at T=4.2K with different 
bias currents 7=0.01-10 mA for (a) the plain and (b) the 
AD bridge. To facilitate the comparison, in (b) the data 
from (a) are shown again as thin lines. 



curve where the slope dR/dH substantially decreases 
(upon increasing i7), can be assigned to the transition 
from the resistive ffux-ffow regime to a resistive regime 
with fully suppressed bulk superconductivity and sur- 
viving ES at Hc2 (and above). This interpretation is 
the same as given in (for surface superconductivity 
with H parallel to the sample surface). However, in 
contrast to our observation, a more gradual transition 
to Rn already at Hc2 without any kinks and no signa- 
ture of ES was observed in 00 when H was applied 
perpendicular to the sample surface. 

We would like to emphasize that the position of the 
kink should be close to the upper critical field Hc2 but 
not identical to it, since the destruction of bulk super- 
conductivity is a thermodynamical property of a ma- 
terial, but the kink can be observed only under strong 
non-equilibrium conditions upon the bias current in- 
jection. Still, below we use the field value where the 
kink appears as the experimentally determined Hc2 
value. 

Obviously, in our case the edge states form contin- 
uous channels with enhanced conductivity, which re- 
duce the overall resistance to a value below Rn- The 
observed reduction of the resistance at the kink fea- 
ture in R{H) with decreasing / can be explained by 
the strengthening of ES upon decreasing /, until at 
small enough currents the injected bias current fiows 
entirely as a dissipationless supercurrent along the 
edge channels at i7 = Hc2^ leading to a disappear- 
ance of the kink feature. 

As described above, the full normal resistance Rn 



is reached (for all / values) at the same field, which 
we now associate with the upper critical field i^cs ^ 
1.7 Hc2 for ES. An analysis of the T dependence of Hc2 
and Hc3 will be presented in Sec. IIII B[ 

In the AD bridge, the holes lead to additional 
"edges" in the sample interior, which results in a 
higher volume fraction of ES and more effective pin- 
ning. This explains the lower R value (as compared to 
the plain bridge) at the kink when bulk superconduc- 
tivity becomes suppressed at i^c2- The origin of the 
second kink at Hc2<\H\<Hc3^ developing at rather 
large bias current [Fig.[2fb)], might be associated with 
a slightly reduced Hcs value at the AD edges, as com- 
pared to the edges of the bridge, due to the different 
edge geometry. However, further investigations are 
required to provide a more conclusive explanation on 
this feature. Similarly, we cannot yet provide an ex- 
planation for the observed shift of the onset of dissi- 
pation to larger iJ, for the AD bridge (as compared 
to the plain bridge) for the lowest values of /. 

To compare the experimental results with theo- 
retical calculations based on the Ginzburg-Landau 
(GL) model described in the Appendix, we calcu- 
lated for a rectangular plain superconducting thin 
film {W = 30^0, L = 60 ^o; is the GL coher- 
ence length at T = 0) the spatial distribution and 
time dependence of the normalized order parameter 
(OP) wave function ?/;(x,?/,t) and the voltage drop 
V{t) along the rectangle for different values of H and 
normalized bias current density j at a reduced tem- 
perature T/Tc=0A7 (corresponds to T=4.2K for Nb 
with Tc=9K). We want to note, that the real dimen- 
sions of the investigated sample exceed considerably 
the dimensions used in our modeling. Nevertheless, 
the model correctly captures the essential physics be- 
hind the discussed effects for H > Hc2- Figure [31(a) 
shows the spatial distribution \ilj{x^y)\ for a rather 
small value of j=5 x 10~^ for five different values of 
H/Hc2 from 0.19 to 1.50. We note that the chosen 
value for j is several orders of magnitude below the 
GL depairing current density jgl =0.386 [cf. the ap- 
pendix] at T=0. In all cases, the OP distributions 
are time- independent ("stationary case") correspond- 
ing to zero resistance. At \H\ < i^c2, a regular vortex 
structure appears and the density of vortices increases 
with increasing H. However, even when bulk super- 
conductivity is depleted at \H\ > Hc2^ superconduct- 
ing channels with finite and time-independent |?/^| run- 
ning along the edges of the rectangle are still present 
and can provide a non-dissipative current transfer. 
If H is further increased, the superconductor turns 
to a non-stationary regime with finite resistance and 
reaches its normal value at the upper critical field for 
ES at = Hc3- The calculated R/Rn vs H/Hc2 
curves for different j are shown in Fig. [31(b). The nu- 
merical simulations reproduce the shift of the curves 
to smaller H and the decrease of the slope dR/dH in 
the interval Hc2 < |i^| < i^c3 as j increases. This is 
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FIG. 3: (Color online) Numerical GL-simulation results 
for a superconducting rectangular thin film {W — 30 ^o, 
L = 60 ^o) biased at normalized current density j at vari- 
able magnetic field H and T/Tc=0.47. (a) Spatial distri- 
bution of the modulus of the normalized order parameter 
wave function |'0(x,^)| with j=0.5xl0~^ [lowest value in 
(b)]. The five panels show simulation results for different 
values of H/Hc2] j is flowing from top to bottom, (b) 
R/Rn vs H/Hc2 for different j. Two vertical dashed lines 
depict the upper critical field Hc2 and the critical field of 
ES Hc2, = 1.695 ifc2. 



qualitatively the same as observed experimentally in 
Fig. [2j However, for large enough j (finite R at Hc2) 
our model is unable to describe the kink in R{H) close 
to Hc2 and the disappearance of Rior H < Hc2 , since 
bulk pinning was not taken into account [cf. curve for 
j = 8 X 10-3 in Fig.[3];b)]). 



B. Superconducting phase diagram for the plain 
bridge 

Figure m^a) shows the results of the R{H) measure- 
ments for T=4.2-8.7K. With increasing T, the devia- 
tion from R = and the kink, both shift to smaller H 
values, and the resistance at the kink shifts to a higher 
R/Rn ratio, while the change in the slope dR/dH at 
the kink becomes less pronounced. In order to ex- 
perimentally determine Hc2 and Hcs we use the field 
value at the kink and a criterion of 0.98 i^n, respec- 
tively. The determined transition lines for Hc2 {T) and 
Hc3{T) for the above mentioned criteria are shown 
in Fig. m^b). The experimental transition line for 
HcsiT) can be fitted with Eq. ([I]) and T^o = 8.5 K 
which extrapolates to i7^3''=15.3kOe. Plotting the 
transition line for Hc2{T) with the relation Hc3{T) = 
1.695 Hc2{T)^ we find that the experimentally deter- 
mined Hc2 values are close to the calculated transition 
line for bulk superconductivity. This result gives con- 
vincing evidence that depending on T, / and H (per- 
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pendicular to the sample surface), our sample can be 
either in the state with developed bulk superconduc- 
tivity and pinned vortex lattice {H < i^c2), or in the 
resistive state, controlled by ES (i^c2 < H < Hcs). 




H (kOe) 



FIG. 4: (Color online) (a) R(H) curves (normalized to 
normal state resistance Rn) of the plain bridge for T=4.2- 
8.7 K (from outside to inside), (b) H — T-phase diagram 
of the plain bridge. Data points for Hcs (dots) and Hc2 
(triangles) are deduced from R{H) curves in (a). The lines 
are the transition lines for ES, which were fitted to the data 
points with Eq. ([T]) and Tco = 8.5 K, and the calculated 
transition lines for Hc2 using the relation Hc3 = 1.695 Hc2 



C. Visualization of the current distribution by 
LTSLM 

We used LTSLM to visualize the current distribu- 
tion in the Nb bridges during the transition from bulk 
superconductivity to the normal state. As the max- 
imum H was limited to ^ 2kOe in this setup, the 
LTSLM measurements were performed at rather high 
rvalues, T = 7.0 -7.5 i^. 

Figure[5](a) shows an i^-series of beam-induced volt- 
age images, AV{x^y), at T=7.5K for various super- 
conducting states of the plain Nb bridge, oriented 
vertically in all these images [cf. optical image (left 
panel) in Fig.[5{a)]. For a more quantitative analysis, 
we show an i^-series of linescans, AV{y), across the 



bridge in Fig. [5jb) and (c). The insets in Fig. [5jb) 

and (c) show the corresponding R{H) curve, from 
which we estimate Hc2 ~l.lkOe and Hcs ^1.8 kOe. 
At H = 0.67 kOe in Fig. [S^a), the LTSLM signal is 
zero, which means that the beam-induced perturba- 
tion is not strong enough to suppress superconductiv- 
ity and to induce a voltage signal. Upon increasing 
the first signal appears at ^ 0.8 kOe which cor- 
responds to the onset of the resistive transition [see 
inset in Fig. E^b)]. With further increasing iJ, the 
signal at the edges is enhanced, but also a signal from 
the inner part of the bridge appears. The latter can be 
attributed to the depletion of bulk superconductivity 
with increasing H (below i^c2), which leads to an in- 
creasing voltage response to the perturbation by the 
laser beam with a maximum beam-induced voltage 
signal at = 1.06 kOe, which is very close to the es- 
timated Hc2 value. The pronounced edge signal below 
Hc2 can be explained by the suppression of the edge 
barrier for vortex entry /exit by the laser spot. Hence 
one can expect that irradiation at the edges of the 
bridge should strongly affect the vortex pattern and 
the resulting current distribution. In contrast, laser ir- 
radiation of the interior of the bridge does not change 
the energy barrier and the modification is probably 
less pronounced and the signal in the interior is much 
smaller. 

Fig. [Sjc) shows linescans for H > Hc2- For fields 
larger than i^c2, the beam-induced voltage in the cen- 
ter of the bridge drops almost to zero while large peaks 
are still observed at the edges of the bridge. This ap- 
parently reflects the fact that above i^c2, the bulk is 
no longer superconducting and therefore does not lead 
to a voltage signal, while the edges still contribute to 
a strong LTSLM signal due to ES. The rather large 
width of these edge peaks in the state of ES can be 
explained by the fact that the edge states are not 
only perturbed when the laser beam spot is centered 
right at the edges, but also when the tail of the beam- 
induced heat distribution leads to a suppression of the 
edge states when the beam is centered slightly off the 
edges. A further increase in H leads to a gradual de- 
crease of the edge peaks which finally disappear at 
H = 1.76 kOe which is close to i^c3- Above i^c3, the 
sample is completely in the normal state and the effect 
of the laser beam on the resistive state is negligible. 

In summary, the linescan series in Fig.[5fc) indicate, 
that above i^c2, the dominant part of the current is 
flowing at the edges of the sample. Thus, LTSLM 
seems to be a capable to visualize the ES states and 
to identify the different regimes in the R{H) curves 
for the plain bridge. 

For comparison, we show a linescan series (variable 
H) for the plain [Fig. [6]; a)] and AD bridge [Fig.lHJb)] 
at T = 7.0 K. The corresponding R{H) curves with 
the bias points of the linescans are shown in Fig.[6l^c). 
The Hc2 value for this temperature is ~ 1.6 kOe. 
We note, that the beam-induced signal of the AD 
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FIG. 5: (Color online) LTSLM signals from plain Nb bridge at T=7.5K. (a) Optical image (left) and voltage images 
AV(x,y) for different H; (b) linescans AV{y) across the bridge for H < Hc2 and (c) for H > Hc2- The insets in (b) and 
(c) show the R{H) curve with corresponding bias points for the linescans. Vertical dashed lines in (b) and (c) indicate 
position of the edges of the bridge 
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FIG. 6: (Color online) LTSLM linescans AV{y) for variable H at T=7.0 K across (a) the plain and (b) the AD bridge; note 
the different AV scale in (a) and (b). Vertical dashed lines indicate position of the edges of the bridge (c) Corresponding 
R(H) curves with bias points for the linescans in (a,b). 



bridge is higher than for the plain bridge, which we 
ascribe to the higher current density in the AD bridge 
due to its reduced cross section because of the holes. 
For the lowest field value, H = 1.20 kOe, the beam- 
induced heating of the laser has no effect, while at 
H = 1.34 kOe the whole cross section of both bridges 
leads to a LTSLM signal. As H increases further, the 
LTSLM signal from the edges becomes larger than the 



signal from the interior and the overall signal increases 
up to = Hc2 = 1.60 kOe. For the plain bridge the 
overall signal gets strongly reduced above Hc2^ and 
the signal from the central part of the bridge almost 
vanishes. The key difference between the plain and 
AD bridge is that for the latter sample the voltage 
signal gets much less reduced and the whole cross sec- 
tion of the AD bridge gives a measurable signal. This 
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means that the current is distributed across the entire 
width of the bridge even for H > Hc2- This observa- 
tion is consistent with the R{H) measurements shown 
in Fig.[6l^c), where the additional edges inside the AD 
bridge lead to a different shape in R{H) and a lower 
R for any value of H within the interval Hc2^H^Hc3. 



within a stationary Ginzburg-Landau model. Since 
the mentioned asymmetry results from the superposi- 
tion of the bias current / and the currents induced by 
the applied H field, the asymmetry can therefore be- 
ing changed either by changing the current direction 
or the sign of H. 



(a) 



D. Bias-current-induced asymmetry: LTSLM 
response and Ginzburg-Landau simulations 

According to Fig.[5fb,c), the LTSLM signal AV{y) 
is asymmetric with respect to the bridge center (axis 
y=0) for several H values around i^c2, i-e. the right 
maximum is slightly higher than the left one. This 
asymmetry in the beam-induced voltage response can 
be explained by an asymmetry in the supercurrent 
density distribution js,x{y) close to the left and right 
edge. Based on the time dependent GL model^, we 
calculate the time-averaged quantities for the OP dis- 
tribution and the x-components of the su- 
perfluid current density {js,x){y) and the normal cur- 
rent density {jn,x){y)' 

Figure [3 shows results of such calculations for H = 
1.3 Hc2 and T/Tco=0.47, which were obtained for zero 
bias current density j = js,x ^ jn,x [Fig. EJa)], for 
j close to the critical current density at = Hc2 
[Fig. [Tl^b)] and for j which is larger than the critical 
current density for ES within the entire field range 
He2<H<H,s [Fig. Etc)]. 

According to our calculations, even in the resistive 
ES state, there is a finite superfluid flow localized 
within the ES channels. These supercurrents are cir- 
culating in opposite direction within each of the two 
edge channels, which is due to the applied magnetic 
field H. 

For further analysis, we determined the net currents 
iLi 'iR and in- Here, il and ir are the integrals of js^x 
across the left and right edge channel, respectively 
(shaded areas in Fig.[7|); in is the integral of jn,x across 
the entire width of the rectangular film. Hence, for 
the normalized bias current = have % = 

ih^in^ in- 

For j = {ih = 0) [cf. FigH^a)], the net currents il 
and iRin the right and left edge channel have the same 
finite amplitude, but differ in sign, and = 0. For 
j > {% > 0) [cf. Fig|7fb,c)] the steady-state distri- 
bution of the superconducting parameters differs from 
the case j = 0. Now, and iR do have the same (pos- 
itive) sign, but different amplitudes. Thus, analyzing 
only the large-scale details in the supercurrent dis- 
tribution (spatially averaged over length scales much 
larger than the coherence length ^o), one can think in 
terms of a combination of two parallel currents flowing 
along the sample edges with different amplitudes de- 
pending both on / and H direction. It should be noted 
that a very similar situation - the asymmetry of the 
critical current density - was described by Park (2Q| 
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FIG. 7: (Color online) Time- averaged normalized OP wave 
function (IV^I^) (2//^o) and x-components of the normal- 
ized superfluid current density {js,x){y/^o) and normalized 
normal current density {jn,x){y /^o), calculated for a rect- 
angular thin film {W = 30 ^o, L = 60 ^o; cf. Fig. [3]) at 
T/Tc=0.47 and H = 1.3Hc2. The three graphs differ in 
the normalized bias current (a) ib=0, (b) ib=0.12 and (c) 
ib— 0.24. iL and ir denote the integrals of {js,x) (shaded 
areas) over the left and right edge channels, respectively. 

In order to prove, whether the asymmetry of the 
LTSLM signal can be related to the bias-current- 
induced asymmetry, we calculated the normalized 
beam- induced voltage Av{y), i.e., linescans across a 
rectangular superconducting thin film with the geom- 
etry as in Fig.[3]and Fig. [71 biased at j = 4x 10~^. De- 
tails of the calculation can be found in the Appendix. 
Assuming a Gaussian shape of the laser-beam-induced 
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increase in T with a maximum amplitude AT and a 
full width half maximum of a = T^o^ we obtain the 
linescan series for different values of H/Hc2 shown in 
Fig. [HI These simulations clearly show that the voltage 
signal has maxima near the left and right edges, and 
that their amplitudes are different, with this asym- 
metry being most pronounced at = i^c2- This is in 
nice agreement with experimental LTSLM results. 
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FIG. 8: (Color online) Calculated normalized LTSLM 
beam-induced voltage IS.v—Vori-Vo^ Ysy/^o across a rectan- 
gular thin film {W = 30^0, L = 60 ^o; cf. Figs. [3] and El) for 
different values of H/Hc2 at T/Tc=0.47 and j=4xl0"^. 
The vertical dashed lines indicate the position of the edges. 




y(Mm) y(|jm) 

FIG. 9: (Color online) LTSLM linescans AV{y) across the 
plain Nb bridge (T=7.2K, ifc2=1.37kOe, |if| = L40kOe, 
|/| = lmA) for different sign of / and (a) negative H and 
(b) positive H. Vertical dashed lines indicate position of 
the edges of the bridge 

The additional edges induced by the holes in the AD 
bridge lead to a different shape of the R{H) curves 
as compared to the plain bridge. The ES state as 
well as the evolution of superconductivity upon sweep- 
ing H was imaged by low-temperature scanning laser 
microscopy (LTSLM). For the ES state, LTSLM re- 
vealed an asymmetry in the currents flowing along the 
left and right edges, depending on the relative direc- 
tion of applied current and external field, as proposed 
long time ago^U. Our calculations based on the time- 
dependent Ginzburg-Landau theory confirm essential 
features of the experimental results. 



To proof experimentally, that the asymmetry de- 
pends on sign of H and /, we performed a series of 
LTSLM linescans on the plain Nb bridge. The rever- 
sal of the asymmetry of the measured LTSLM signal 
upon the inversion of the / and H signs is illustrated 
in Fig. ini^a) and (b). We find that the right peak is 
larger for / > while the left peak is larger for / < 
and vice versa. The slightly larger amplitudes of the 
peaks in Fig. [9l^b) are probably due to the residual 
field (in the 10 Oe range) in the cryostat at the sam- 
ple position. To the best of our knowledge, this is the 
first direct experimental verification of an asymmetry 
in the current density in the ES state, as predicted by 
Park for surface superconductivity. 



IV. CONCLUSION 
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In this paper we studied experimentally and nu- 
merically the peculiarities of the resistive transition 
in thin-film Nb microbridges with and without anti- 
dots (ADs) in perpendicular magnetic field H. From 
integral R{H) measurements we find that the transi- 
tion from bulk to edge superconductivity (ES), and 
finally to the full normal state, can be identified by 
a pronounced change in slope dR/dH^ which however 
strongly depends on the applied bias current density. 



VI. APPENDIX 

In order to describe the general properties of the re- 
sistive state in a mesoscopic superconducting thin film 
sample and to compare them with experiment, we use 
a simple time-dependent Ginzburg-Landau (TDGL) 
model. For simplicity we assume that the effect of 
the superfluid currents on the magnetic field distribu- 
tion is negligible and consider the internal magnetic 
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field B equal to the external magnetic field H (per- 
pendicular to the thin film plane). This assumptions 
seems to be valid for the following two cases; (i) for 
mesoscopic thin-film superconductors with lateral di- 
mensions smaller than the effective magnetic penetra- 
tion depth A = X^/d (Al is the London penetration 
depth, d is the thickness); (ii) for superconductors for 
large H and/or T (i.e. close to the phase transition 
line), when the superfluid density tends to zero. Then 
the TDGL equations take the form 

u(^^^^i^^i^ = r {i^- \^\^^) + (V + iAf (2) 
r = l-T(r)/T,o, (3) 
VV = div j., js = -\t {v^* (V + zA) - c.c.}, (4) 

where i\) is the normalized order parameter (OP), ^ 
is the dimensionless electrical potential, A is the vec- 
tor potential [rot A = i^e^], T(r) is local tempera- 
ture (potentially position-dependent), js is the den- 
sity of the supercurrent, u is the rate of the OP re- 
laxation, c.c. stands for complex conjugate. We use 
the following units: m*crn/3/(2e^a) for time, the co- 
herence length ^0 at temperature T = for distances, 
$o/(27r^o) for the vector potential, fie|a|/(m*crn/3) for 
the electrical potential, and 4ea^^o/(^/^) for the cur- 
rent density, where a = — ar and /3 are the conven- 
tional parameters of the GL expansion, e and m* are 
charge and the effective mass of carriers, Gn is the nor- 
mal state conductivity. In these units the Ginzburg- 
Landau deparing current density at T = is equal 
to 0.386. We apply the boundary conditions in the 
following form 



d_ 
dn 



^lAn] = 0, 



dip 
dn 



■ Jext-) 



(5) 



where n is the normal vector to the sample's bound- 
ary F, jex^ is the normal component of the inward 
(outward) flow of the bias current density j (with 
\j\ = j). We do not consider bulk pinning, since the 
number of additionally required parameters (describ- 
ing the spatial distribution of pinning sites and their 
pinning strength) would be too large. 

We calculate^^ the instant value of the normalized 
voltage drop v{t) = {(pi{t)) — (9^2 (^)) and analyze the 
dependence of v{t) on H and jext- Here 



ipi{x,y,t)dxdy 



(6) 



is the time-dependent electrical potential averaged 
over the region Si ("virtual electrodes", i = {1,2}). 
These regions have the same width as the sample 
width and they are shifted from the physical edges 
towards the sample interior (see Fig. [TQ|) for eliminat- 
ing the effect of the sample edges. In addition we for- 
mally consider an inhomogeneous sample, containing 



two areas at the left and right edges [cf. Fig.[TO] with a 
critical temperature Td (at H = 0) and upper critical 
field H^2 1 (at T = 0) exceeding considerably Tco and 
H^2 ^^^^ of the sample. The reason for that is 

a pure technical one. This approach guarantees that 
the injected normal current is fully converted into a 
supercurrent within these enhanced superconducting 
areas at any temperature and any value of H. 



^ Regions for the calculation of 
the electrical potentials 

<(p,(t)> and <(p2(t)> 
averaged over S, and S2 



FIG. 10: (Color online) Schematic drawing (top view) of 
rectangular superconducting thin film bridge considered 
for GL simulations. Arrows indicate injection and extrac- 
tion of the bias current. Shaded areas Si and S2 are virtual 
electrodes. 



For the stationary regime all the calculated param- 
eters, after transient processes induced by changes 
in the external parameters, tend to their time- 
independent values, pointing out to the absence of en- 
ergy dissipation for the established state and R ^ 0. 
For larger T, H or the relaxation to the station- 
ary case becomes impossible and all parameters oscil- 
late in time. Calculating the mean normalized voltage 
drop V, averaged over a very large time interval (in- 
cluding up to 10^ of the voltage oscillations), one can 
determine the normalized beam-induced LTSLM volt- 
age signal A'^ = — Vq^^ , where Von and '^off are the 
time averaged normalized voltage signals if the laser 
beam is on or off, respectively. 

The effect of the focused laser beam can be treated 
as a quasistatic perturbation of the superconducting 
properties of the bridge, since the time scales of this 
perturbation are much longer than the GL time con- 
stant. In the most simple form this perturbation can 
be modelled as a Gaussian-like increase in local tem- 
perature in Eq. (j3j): 



T(r) = To + AT • e 



[-ix-xof-iy-yof]/cT^ 



(7) 



Here, Tq is the sample temperature if the laser beam is 
off or far from the beam spot centered at (xq, ^o), AT 
is the amplitude of the local heating, depending on 
the beam intensity and on the rate of heat dissipation 
due to the thermal conductivity of the superconduct- 
ing film and the substrate and on the thermal bound- 
ary resistance between the film and the substrate; a 
is the full width half maximum of the beam-induced 
temperature profilei^ 
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